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We study the collision probability p of particles advected by open flows displaying chaotic advec- 
tion. We show that p scales with the particle size S as a power law whose coefficient is determined 
by the fractal dimensions of the invariant sets defined by the advection dynamics. We also argue 
that this same scaling also holds for the reaction rate of active particles in the low-density regime. 
These analytical results are compared to numerical simulations, and we find very good agreement 
with the theoretical predictions. 
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Many fluid flows of interest to science and to engi- 
neering are open flows, which are characterised by the 
presence of inflow and outflow regions [l| (see Fig. [TJi). 
The dynamics of advected particles in open flows is tran- 
sient, with typical particles leaving to the outflow region 
in finite time. Advection in many open flows is chaotic 
@; 13 ; examples of chaotic open flows are found in areas 
such as microfluidics climatology 0L physiology @, 
population biology Q and industry The dynam- 

ics of advected particles is governed by the chaotic sad- 
dle which is a set of unstable orbits contained in a 
bounded region of space known as the mixing region (see 
Fig. [T^,). Fluid elements are repeatedly stretched and 
folded in the mixing region, and this causes any portion 
of the fluid to be deformed by the flow into a complex 
filamentary shape, which shadows the unstable manifold 
of the chaotic saddle, defined as the fractal set of orbits 
which converge asymptotically to the chaotic saddle for 
t -» -oo fioj. 

The fractal structures generated by chaotic advection 
have dramatic consequences for the dynamics of active 
processes taking place in open flows, such as chemical re- 
actions and biological processes The stretching and 
folding of fluid elements by the flow tends to increase 
the area (or perimeter, in 2D flows) of contact between 
two reacting species, which results in an enhancement of 
the reaction due purely to the advective dynamics of the 
flow — this has been called dynamic catalysis Q. It has 
been shown [j| [ll[ that this effect appears as a singular 
production term in an effective reaction rate equation, 
which has a power-law dependence on the amount of re- 
actant; the coefficient of the power law is determined by 
the information fractal dimension of the unstable mani- 
fold of the chaotic saddle. Later works have established 
that this enhancement of activity by chaos is a very gen- 
eral and robust phenomenon Q , and is a feature of many 
kinds of activity and flows, including non-periodic [12j, 
non-hyperbolic [l3, 14 1 and three-dimensional [l5[ flows. 

All the works so far on the effects of chaos on chem- 
ical or biological activity in open flows have adopted 
continuum descriptions of the reacting material — ei- 
ther describing the reactants using continuous concen- 



tration fields, or modelling the propagation of the reac- 
tion through reaction fronts. This approach ignores the 
fact that any activity taking place in the fluid ultimately 
arises from collisions of particles. If the number of react- 
ing particles is very large, so that the mean free path be- 
tween collisions is much smaller than the typical length 
scale of the system, then this continuous description is 
expected to be valid. However, in the low-density limit, 
where typically particles traverse large distances before 
colliding and reacting, the concepts of reactant concen- 
tration and reaction front have no meaning. In this case, 
one must adopt a kinetic theory approach, where activ- 
ity is described in terms of the probability of collision 
between particles. Although this idea has been used in 
closed flows — in the context of chemical catalysis [l6| , 
particle coalescence [l7 , 18 1 and crystallisation [nj], for 
example — so far this idea has not been pursued in the 
case of open flows, despite the importance of the low- 
density regime for applications, which is appropriate for 
describing systems as diverse as platelet activation in 
blood flows plankton population dynamics (l^ . and 
raindrop formation 2^, 21 1. 



In this work we investigate the effect of chaotic ad- 
vection on the activity of particles in the low-density 
limit. The rate of reaction in this limit is determined 
by the collision probability p(S) of two particles com- 
ing within a distance d of each other before they escape 
to the outflow; we assume a reaction event takes place 
when that happens. 5 can be thought of as the size of 
the particles, or alternatively as determined by the re- 
action cross-section a (with 6 ~ \f5)- We show that p 
scales with 8 as a power law, p(S) ~ 6? , and we derive an 
analytical expression for the coefficient (3 in terms of the 
fractal dimensions of the stable and unstable manifolds 
of the chaotic saddle. The classical kinetic theory ap- 
proach, based on the Smoluchowski equation (22| . which 
assumes there is a homogeneous mixing of particles in 
the mixing region, would predict j3 = N, where N is 
the spatial dimension, but we show that because most 
collisions happen in the neighbourhood of a fractal set, 
Smoluchowski's result does not hold, and (3 ^ N for open 
chaotic flows. In the low-density limit, the total number 
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FIG. 1: (a) Illustration of an open flow, (b) Schematic 
depiction of the stable (W s ) and unstable (W u ) 
manifolds in an open flow. The grey bands are 
8- neighbourhoods of W s and W u (see text). R\ and R2 
indicate regions where initial conditions are chosen. 



Q{8) of reaction events between reacting particle species 
which takes place in a given flow also scales as 8" . We 
argue that Q is a good measure of the efficiency of mixing 
in open flows, and thus (3 characterises the mixing of an 
open flow. These analytical predictions are compared to 
numerical simulations in 2D flows, and we find the results 
agree very well with our theory. 

We start by deriving the scaling of the collision prob- 
ability p{8) of two particles with randomly chosen initial 
conditions, by first assuming that after the particles en- 
ter the mixing region (Fig. [TJt), the flow acts as a perfect, 
uniform mixer, until the particles escape to the outflow. 
This perfect mixing assumption means that the positions 
of the two particles are effectively randomised in a time 
interval r which is characteristic of the flow. The proba- 
bility that two particles randomly placed in a bounded N- 
dimensional region are within a distance 8 of each other 
scales as 8 , Assuming that 8 is small enough so that 
the probability p(S) that the two particles collide before 
escaping satisfies p(8) <?C I, we conclude that p(8) is pro- 
portional to S N and to the average time T the particle 
stays in the mixing region before escaping. p(S) therefore 
scales with 8 as 



p(8) ~ 8 



(1) 



The perfect mixing assumption is the assumption used 
in Smoluchowski's classical coagulation theory (22[, and 



we will refer to this 8 N scaling as the Smoluchowski pre- 
diction. 

The Smoluchowski scaling derived above does not take 
into account the fractal nature of the particle distribu- 
tions generated by chaotic advection in open flows. We 
will now derive the correct scaling of p{8). Let two par- 
ticles be chosen randomly within two bounded and dis- 
joint regions R\ and R2 (see Fig. [TJd). As a result of the 
chaotic dynamics of the flow in the mixing region, those 
trajectories which take longer to escape converge towards 
the unstable manifold W u . Consequently, the particles 
will only have a. chance to meet if both come within a 
distance 8 (the reaction distance) of the unstable man- 
ifold. This will take place if their initial conditions are 
within a distance of order 8 of the stable manifold W s , 
which is the set of initial conditions converging to the 
chaotic saddle for t — > 00 (see Fig. \Tjp) The proba- 
bility pP(S) that one particle with initial conditions in 
i?i lies within a distance of order 8 of the stable manifold 
is proportional to the volume V± (or area in two dimen- 
sions) of the intersection of Ri and a (^-covering of the 
component of W s where the conditionally invariant mea- 
sure is concentrated fiol . 23 1 (see Fig. QJ>). If one imagines 



i?i to be covered by a grid of size 8, the number of cells 
in the grid intersecting the component of W s where the 
conditionally invariant measure is concentrated scales as 
N x {8) ~ <5- D i [13, where Df is the information dimen- 
sion of the stable manifold. Since the volume of each cell 
in the grid is 8 N , the volume of the portion of the grid 
which intersects W s scales as V Y ~ N(8)8 N ~ S n ^ d k 



We therefore have Pi ; (8) 
ing holds for the probability P< 

is close to W s . Therefore, the probability P s = P^P^ 
that both particles reach a (^-neighbourhood of the un- 
stable manifold scales as 



S N-D 
(2 



and the same seal- 
that the other particle 



Ps(S) 



c2N-2D? 



(2) 



The two particles collide if they reach the unstable mani- 
fold within a distance 8 of each other. Considering again 
a grid of size 8 covering the unstable manifold W u , the 
probability P u (8) that two particles collide can be ap- 
proximated by the probability that they end up in the 
same cell in the grid. This is given by P u {8) = X^M?, 
where pi is the (conditionally invariant) measure of cell 
i, and the sum is over all cells intersecting the unstable 
manifold. From the definition of the correlation dimen- 
Tol . 23 1, we find that in the limit 8 — > we get 



sion D% 



Pu{8) 



(3) 



The overall probability p(S) — P S (8)P U (S) that the two 
particles will collide therefore scales as 



p(S)~8 , with (3 = 2N-2Dl 



(4) 



This is our main result. This scaling is clearly different 
from the Smoluchowski prediction: for chaotic flows, in 




FIG. 2: (a) Numerical approximation of the unstable 
manifold for the blinking vortex-sink flow, with rj = 0.5 
and £ = 10; for these parameters, the dynamics is 
chaotic and hyperbolic, (b) Location of collisions for 
S = 10 -3 , for the flow with the same parameters as 
above. It is clear that collisions take place in the 
vicinity of the unstable manifold. The initial conditions 

were chosen in the regions R\ (—0.5 < x < 0.5, 
< y < 0.15) and R 2 (-0.5 < x < 0.5, 0.2 < y < 0.35). 
The time step used in the simulations was At = 0.02. 



general we expect j3 ^ N . However, (3 does approach the 
Smoluchowski value of N in the limit where the escape 
rate vanishes, when D 2 ,D 2 — > N, and from Eq. (0| this 
implies that /? — > N. This behaviour is expected, since 
in this limit the dynamics within the mixing region can 
be considered as that of a closed container being stirred, 
with a small leak, when particles spend a long lime in 
the mixing region and the perfect mixing assumption is 
approximately valid. 

The scaling predicted by Eq. @ is valid for any choice 
of regions Ri and R 2 , as long as both regions intersect 
the stable manifold W s . If there is no intersection, the 
two particles never meet because they do not get to the 
mixing region, and collisions cannot take place. 

To test the prediction of Eq. (H}, we will use the blink- 
ing vortex-sink flow, [zj . |25| which is a generalisation of 



FIG. 3: Collision probability (squares) and the number 
of reaction events (circles) as a function of the reaction 
distance S for the blinking vortex-sink flow, with the 
same parameters as in Fig. [5] The initial conditions 
were chosen in the regions R\ and R 2 defined in Fig. [21 
The number of pairs used varied with 6, and it was such 
that it guaranteed that at least 10 3 collisions took place 
for each 5; typically this meant the number of pairs was 
of order 10 10 -10 11 . The straight lines are the fitting of 
the data to a power law of the form 5@ ; in both cases 
fitting yields j3* = 2.23, compared to the prediction 
P = 2.25 of Eq. dU). The collision probability for the 
area-preserving baker map with escape is also shown 
(diamonds), for contraction factors 0.3 and 0.45 27|. 



Fitting yields /3* = 2.30. For this system Df = 1.701 
and D 2 = 1.698, and from Eq. Q the prediction is 
/3 — 2.298. The triangles show the collision probability 
for the conservative Henon map for the nonlinearity 

parameter A = 7 for which Df = 1.58 and 
D 2 = 1.32, leading to a prediction of /? = 2.16; fitting 
yields (3* — 2.14. The baker map data (diamonds) in 
this plot has been shifted in the graph above to improve 
visibility. 



Aref's blinking vortex flow [26]. It is a 2D incompress- 
ible flow on an infinite plane, with two sinks at positions 
(+a,0) and (— a,0), which open and close periodically 
in alternation: in the first half of each cycle one sink is 
open and the other one is closed, and in the second half 
the situation is reversed. Each vortex-sink is modelled 
as a point source of vorticity superimposed to a localised 
sink, and flow which falls on either of the sinks disappears 
from the system and does not come back. This system is 
clearly an open flow, where the inflow region corresponds 
to the whole space beyond the sinks. 

Assuming that the particles advected by the flow be- 
have as passive tracers, their equations of motion are [24j 



r = —C/r, if = K/r 2 



(5) 



where r and ip are polar coordinates whose origin alter- 
nates between (a, 0), during the first half of each period, 
and (— a, 0), during the second half of each period; C 
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and K are parameters representing the strengths of the 
sink and the vortex, respectively. The advection dynam- 
ics is determined by the two dimensionless parameters 
77 = CT/a 2 and £ = K/C, where T is the flow period. 
We choose units so that a = 1 and T = 1. The dynam- 
ics is chaotic for a wide range of 77 and £ (24j. Fig. [^Ja. 
shows the unstable manifold for one particular choice of 
parameters, 77 = 0.5 and £ = 10. 

We computed the collision probability p(S) by follow- 
ing the trajectories of many pairs of particles until one 
particle in each pair escapes, and at every time step At 
we check if the distance between the two particles is less 
than 8; if so, the particles are considered to have collided. 
The initial condition of one of the particles in each pair is 
chosen randomly within region Ri , defined in the caption 
of Fig. [2j and the other particle in the pair is started in 
region R2 (see Fig. [5]). The fraction of the pairs which 
collide before escaping gives us a numerical approxima- 
tion of p(S). Figure [3] (the squares) shows the result of 
applying this procedure for a range of values of 8. The 
resulting scaling is clearly a power law p(S) ~ 8@ , and 
fitting yields j3* ps 2.25. To compare this with the predic- 
tion of Eq. ((3]), we use the value D\ ps 1.74 from [24]] ; to 
compute D2 we first find an approximation W*? 1 of the 
unstable manifold using the sprinkler method [10|. and 
then we compute the scaling with e of the total number 
of pairs of points on W^ pr separated by a distance less 
than e; the coefficient of the resulting power law is 
[23l | . We found D 2 ps 1.71. This is a two-dimensional 
system, and so N = 2 in Eq. (g]). Equation (g]) then 
predicts j3 = 2.23, which agrees with the value j3* = 2.25 
from the simulation, to within numerical errors. 

Fig. [5^ is a snapshot of W u at the start of a period; 
the shape of W u changes periodically with time, with the 
same period as the flow. Figure^ shows the positions of 
the collisions that take place at the beginning of each pe- 
riod: each dot in Fig. [^b is one such collision. Comparing 
Figs. [2Ja. and^, it is clear that the collisions take place 
in a small neighbourhood of the unstable manifold W u , 
confirming that the assumption made in the derivation 
of Eq. (g]) is correct. 

We have analysed other choices of the parameters 77 
and £ for which the flow is chaotic, and we always find 
that the prediction from Eq. (g|) and the results from 
the simulations differ by less than 1% in all cases. We 
found that the choice of the time step At does not affect 
the scaling of p(8) with 8 — although the value of p for 
any given 5 of course decreases with At, since with larger 
time steps it is more likely collisions will be missed by the 
simulation. 

We have also studied the dynamics of collisions in other 
two-dimensional dynamical systems, including a variant 
of the conservative Henon map which displays transient 
chaos and a version of the baker's map with escape 
[27l | . In all cases, and for all parameter choices, we found 
excellent agreement between the numerical value of /3 and 



the prediction of Eq. (g) (see Fig. 0). 

We now examine a system of reacting particles ad- 
vected by a chaotic open flow. Taking as an illus- 
trative example the case of an acid-base-type reaction 
A + B — > C, consider particles of two species, which we 
will label simply species A and B, such that when an A 
particle comes within a distance 8 of a particle of type 
B, a reaction occurs, resulting in a particle of type C, 
which is the product of the reaction. If we throw Ni 
type- A particles in a region R\ of the flow, and N 2 type- 
B particles in region R2, then we expect that the total 
number N p of type-C particles produced in the system is 
proportional to NiN2p(5) — as long as N\ and N2 are 
not too large, and the low-density assumption holds. We 
therefore predict that the amount N p of product gener- 
ated by the reaction scales with the reaction distance 8 
as with p given by Eq. (gj. 

We put this prediction to the test in the blinking 
vortex-sink system, using an efficient algorithm for find- 
ing the neighbours within a distance 8 [28j . The resulting 
scaling of N p with 8 is shown in Fig. [3] We find a power- 
law scaling as expected, and fitting yields the coefficient 
f3* = 2.25, which is again within 1% of the predicted 
value of j3 = 2.23. We found the same agreement for 
other choices of parameters for the blinking vortex-sink 
system, and also for the other dynamical systems we in- 
vestigated. 

We note that, even though we used the reaction A + 
B — > C as an example above, all kinds of reactions and 
active processes involving particles will be affected by the 
scaling (g|), since all reactions require particles to come 
close together, and the probability of that happening is 
governed by the coefficient /3 of Eq. (gj . We propose that 
the collision probability p, and the way is scales with the 
collision distance 8, is the natural quantity to characterise 
the efficiency of mixing of an open flow (an alternative 
definition is given in [9|). Equation (g]) shows that the 
mixing efficiency is determined by the fractal properties 
of the invariant sets associated with the chaotic saddle — 
namely its stable and unstable manifolds. The discussion 
above shows that the mixing efficiency in turn determines 
the efficiency of chemical reactions (or, in general, active 
processes) in open flows, which means that the "chemical 
efficiency" of a flow is also measured by the coefficient (3 
defined in Eq. (g]). 
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